Localization properties of groups of eigenstates in chaotic systems 
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In this paper we study in detail the localized wave functions defined in Phys. Rev. Lett. 76, 
1613 (1994), in connection with the scarring effect of unstable periodic orbits in highly chaotic 
Hamiltonian system. These functions appear highly localized not only along periodic orbits but also 
on the associated manifolds. Moreover, they show in phase space the hyperbolic structure in the 
vicinity of the orbit, something which translates in configuration space into the structure induced by 
the corresponding self-focal points. On the other hand, the quantum dynamics of these functions 
are also studied. Our results indicate that the probability density first evolves along the unstable 
manifold emanating from the periodic orbit, and localizes temporarily afterwards on only a few, short 
related periodic orbits. We believe that this type of studies can provide some keys to disentangle 
the complexity associated to the quantum mechanics of these kind of systems, which permits the 
construction of a simple explanation in terms of the dynamics of a few classical structures. 
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I. INTRODUCTION 

The investigation of the quantum manifestations of 
classical chaos is at present a very active field of research 
0. The relationship between quantum dynamics and 
classical invariants, mainly periodic orbits (PO) and their 
associated manifolds, is not satisfactorily understood yet 
In this respect, an important achivement would be 
to disentangle the time evolution of a quantum system 
whose classical analog is chaotic, explaining this dynam- 
ics in terms of classical structures. 

Fifteen years ago, a seminal step was done by Heller @, 
who considered the evolution of a wave packet localized 
on an unstable PO. He found that the dynamics in the 
short time limit are controlled by the motion around the 
orbit and its Lyapunov exponent and is well understood 
in terms of semiclassical expressions . 

For longer times the dynamical mixing induced by the 
chaotic behavior of the motion makes the situation more 
complicated, and one single orbit and the linearized dy- 
namics around it is not enough to account for the evo- 
lution of the wave packet. Some time later, Tomsovic 
and Heller || showed that even for those long times, for 
which classical hne structure had developed on a scale 
much smaller than h, the semiclassical propagation of the 
packet can be carried out with good precission, by com- 
puting the corresponding correlation function, C sc i(t), as 
a sum of contributions of the homoclinic excursions of the 
PO. This procedure has nevertheless some drawbacks, 
since for example the number of orbits that is necessary 
to include in the calculation to obtain converged results 
grows dramatically with time. 

A different approach, that maintains an interpretation 
of wave mechanics based on the more simple classical ob- 
jects represented by POs has been presented elsewhere 
In this theory || all quantum information of a 



bounded chaotic Hamiltonian system can be obtained us- 
ing a number of POs that grows only linearly with the 
Heisenberg time. 

Another point worth considering is the importance 
and/or convenience of using wave functions which are 
"dynamically adapted" in this type of studies . By "dy- 
namically adapted" (8j we mean structures that contain 
information not only about the PO, but also about the 
linearized dynamics around the orbit. When considered 
in phase space, these structures appear no only localized 
over the PO, but are also spread over the corresponding 
manifolds. 

The construction of functions of this kind can be car- 
ried out using different approaches. In Ref. |j) the neces- 
sary information is obtained by Fourier transforming the 
short-time exact quantum dynamics, so that wave func- 
tions averaged over the PO path are obtained. Kaplan 
and Heller [|l0| used a sum of coherent Gaussian wave 
packets centered along the PO to obtain the desired func- 
tion. Finally, Vergini and Carlo [^1) constructed semi- 
classically resonances along POs with minimum energy 
dispersion. 

In spite of the differences existing among these func- 
tions, they are all suitable to efficiently investigate the 
quantum manifestations of the classical phase space com- 
plexity in classically chaotic systems. 

The purpose of this paper is threefold. In the first 
place, we focus on the analysis of the topological char- 
acteristics (i.e. distribution of quantum probability) in 
phase space of this type of wave functions highly local- 
ized on an unstable PO. In the second place, special at- 
tention will be paid to the relationship existing between 
the quantum dynamics of this function and that of other 
classical POs of the system located in their vicinity, in 
order to explore the existence of dynamical connections 
and active interplay between them. Finally, another Ln- 
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teresting aspect of these functions is their revivals, that 
take place in some cases for a surprisingly large number 
of cycles without showing any visible sign of dispersion 
|[I2|| . This fact, which is particularly surprising in view of 
the highly chaotic character of the problem that will be 
considered, is due to the special dynamical properties of 
these functions. Summarizing, it should be stressed that 
the most important result of our work is that, proceed- 
ing in this way, only a few short POs are necessary to 
adequately explain the quantum dynamics of this type of 
systems. 

For our study we will use localized functions as defined 
in Ref. J9|, applying them to the study of the dynamics 
in the Bunimovich stadium billiard |l3| , a system which 
is known to behave ergodically from the classical point 
of view. 

The organization of the paper is as follows. In Sect. |l| 
we describe in detail how the method of Ref. || can 
be used to construct the desired wave functions, highly 
localized on a given unstable PO and their associated 
manifolds. Numerical results concerning the topological 
characteristics, both in configuration and phase space, of 
these functions and the corresponding time evolution will 



be presented and analyzed in Sects. HI and IV, respec- 
tively. Finally, the main conclusions derived from our 
work are presented in Sect. 0. 



II. CONSTRUCTING LOCALIZED WAVE 
FUNCTIONS 

In this section we briefly describe the method Q that 
will be used to construct states initially located on a given 
PO and the corresponding manifolds. The procedure has 
been previously applied to the analysis of the scarring ef- 
fect of unstable POs in chaotic systems [ p4[ , as well as 
to the study of the quantum dynamics that takes place 
beyond the first recurrence of the corresponding motion 
0. As it will be shown, these initial states can be used 
as an important tool to disentangle the complexity of the 
full quantum dynamics taking place in the area of influ- 
ence of a given PO, making use of the associated classical 
mechanics as a guide in the process. 

Initially, the method considers a wave function, \<fi(0)), 
well localized in phase space in the vicinity of a given 
PO. The corresponding time evolution can be followed 
by means of the associated autocorrelation function, 



C(t) = (<P(0)\m = <0(O)|e^ 4 |0(O)). 



(1) 



Recurrences in C(t) are known to determine the low res- 
olution structure of the corresponding spectra [Q , 
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and are the origin of the scarring effect, as first discussed 
in ||. Wr(t) is a window function to take into account 
possible different resolutions in It(E). 



A widely spread choice for 0(0)) is a minimun uncer- 
tainty harmonic oscillator coherent state (h is set equal 
to one throughout this paper) 
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exp [tiffo - r°)] , 



(3) 



where (r°,P°) represents the coordinates and conjugate 
momenta of a suitable phase space point along the se- 
lected PO. 

For bounded systems, the evolution of this packet can 
be followed quite conveniently by projection on a com- 
plete set of eigenf unctions, |n), of the Hamiltonian, H , 
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where {E n } is the corresponding set of eigenvalues. In 
our case the billiard eigenstates have been computed us- 
ing the scaling method [^5) . 

Regarding the smoothing function Wt (t) one has many 
choices; i.e. "hat" , exponential or Gaussian functions, 
which render sinx/x, Lorentzian or Gaussian-type line 
envelopes, respectively. In our case we use the Gaussian 
decay: 
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for which the spectrum takes the form 
1 
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(5) 



(6) 



The center of the initial wave packet defined in Eq. (g) 
will follow for some time the classical PO on which it was 
launched, with the dispersion rate evolving according to 
the Lyapunov exponent of the orbit [Q. 

In a second step B , a wave function highly scarred by 
the PO can be constructed very efficiently by consider- 
ing the average of \4>(t)) along the dynamics generated 
by the PO during a given amount of time, usually taken 
of the order of the PO period. Let us remark the im- 
portance of this last requirement, since if the packet is 
allowed to evolve for a very long time, so that due to the 
chaotic character of the dynamics it has the opportunity 
to sample all the available phase space, just the (compli- 
cated) eigenfunctions of the system will be obtained. The 
corresponding expression for the localized wave function 
corresponding to an energy Eq (usually taken as the cen- 
ter of one of the bands appearing in the low resolution 
spectrum generated by \<fi(0))) reads as follows: 
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Phase space representations for the wave functions de- 
scribed in this section can be constructed in a number 
of different ways. In this paper we will follow the pro- 
cedure described in Ref. |lq| , which relies on the use of 
normal derivatives of the wave functions evaluated at the 
boundary of the billiard. On this boundary, Birkhoff co- 
ordinates J!?]] , (q 7 p) , are used to define both classical and 
quantal Poincare surfaces of section, such that g, is the 
arc length coordinate, and p = p • t/|p| is the fraction of 
tangential momentum. The coherent states necessary to 
construct this representation (Husimi functions |l8|]) are 
then defined as: 
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x exp [ip°(q — q )]. 



(8) 



This expression corresponds to a boundary wave packet 
at point (q°,p°) in the surface of seccion; representing a 
bounce off a given boundary point with a specific tangen- 
tial momentum. Then, for a wavefunction with normal 
derivative on the boundary ?/>(<?) (extended periodically 
to the real line), the corresponding Husimi function is 
given by 



H(q°,p°) = \(G q o p0 \iP) 



(9) 



III. TOPOLOGICAL CHARACTERISTICS OF 
THE LOCALIZED WAVE FUNCTIONS 

The system that we have chosen to study is a particle 
of mass 1/2 moving in a desymmctrizcd stadium billiard 
with Newman boundary conditions on the symmetry axis 
(only even-even parity wave functions will then be con- 
sidered). The radius r is taken equal to unity and the 
enclosed area is l+ir/4. This system is known to exhibit 
a great degree of chaoticity, both from the classical and 
quantal points of view. 

In this case, and due to the symmetry of the problem, 
symmetry adapted initial wave functions should be used; 
the corresponding expression being given by 

(x,y\(f>(0)) = G x o^ y opo P a + G_ x a y a_pa P o 



G , 



c.c. 
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is obtained by imposing Newman boundary condition at 
the symmetry axis. 

We consider the dynamics dominated by the horizon- 
tal PO, running along the x axis with y = 0. Accord- 
ingly, we use for our study a symmetry adapted Gaus- 
sian wave packet [Eq. (foj)] with centers defined by the 



phase space point (x , y 



p o pen 



(1,0,A:,0) and width 



(Ti; = (Ty = 1.603/fc 1 / 2 , where k is the usual wave number. 
The corresponding autocorrelation function for a packet 



with an energy at the center of £'=3600, for which the pe- 
riod of the horizontal PO is equal to Xo=l/30, is shown in 
Fig. |l|. As can be seen, two main peaks exist at t ~ 0.016 
and 0.033 respectively, being the first one substantially 
lower than the second. After that, the correlation func- 
tion becomes very complicated. The relative intensity of 
the two peaks can be easily explained if one takes into 
account that the initial coherent state that we are using 
is the sum of two packets, launched with opposite values 
of the momentum. The two packets first separate, one 
being scattered by the vertical axis and the other by the 
circular hard wall of the billiard, so that when they first 
return to the initial point, at a time roughly equal to 
To/2, they have picked up different phases. This gives 
rise to the first (small) recurrence peak in \C(t)\. After- 
wards, the two components continue moving, colliding 
with the circular and vertical walls, respectively. In this 
way, when some time later, at t ~ To, they return again 
close to the initial point, both arrive with approximately 
the same collisional phase, this giving rise to the second, 
more pronounced peak in the autocorrelation function. 
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FIG. 1. Modulus of the autocorrelation function corre- 
sponding to a symmetrized wave packet [Eq. (|^)] initially cen- 
tered on the horizontal periodic orbit of the desymmetrized 
stadium billiard with Newman boundary conditions on the 
axis at an energy _E=3600. The radius is taken equal to 1, 
the enclosed area l+7r/4, and the mass of the particle 1/2. 

Let us consider next the localization of the wave func- 
tions obtained with the method described in the previous 
section [see Eq. (0)] for different values of the smoothing 
time. We choose T=0.01, 0.03, 0.06, and 0.1, in order to 
cover a range of values smaller and greater than To. The 
results are shown in Fig. where we present the cor- 
responding squared wave functions, Husimi based quan- 
tum sufaces of sections, and the associated portions of 
the spectra generated from |0(O)}. Three comments are 
in order. 
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FIG. 2. (Top) Localized (squared) wave functions calculated using Eq. (|7j), (middle) Husimi based quantum Poincare surface 
of section, and (bottom) spectra (Too as a stick spectrum and It in solid line) for different values of the smoothing time, T: 
(a) 0.01, (b) 0.03, (c) 0.06, and (c) 0.1 at _Eo=3700, corresponding to the autocorrelation function of Fig. jjj. The invariant 
manifolds of the horizontal periodic orbit have also been drawn superimposed in the middle pannel. 



In the first place, the wave functions in configuration 
space (top tier) appear localized along the PO, present- 
ing a number of nodes, m, in agrement with the value 
obtained from the Bohr-Sommerfeld quantization condi- 
tion 

k= T (rn+~) (11) 

when applied to the spectral peaks (bottom tier). For 
the horizontal PO, L = 4 and the Maslov index z/=3. 

In the second place, and much more importantly, in 
the localized wave functions it is clearly observable the 
distorsion in the structure induced by the self focal point. 
A discussion of this effect can be found in Ref. [0 . This is 
a very important point, and is a consequence of the fact 
that these functions have been defined in such a way- 
averaging over the motion along the orbit-that they are 
not only located in the right position of configuration 
space, but contain also information on the dynamics of 
the system in the vicinity of the PO. This is more clearly 
seen when examining the phase space picture of these 
functions, i.e. middle tier of Fig. g. There, it is observed 
that the probability density do not just localize over the 
fixed point corresponding to the scarring orbit, but it is 
also spread significantly along their manifolds. 



Finally, this structure over the manifolds changes as 
the value of T is increased, giving the original packet the 
opportunity to explore more and more of the dynamics 
induced by the PO. Indeed, we see in Fig. || that for 
T=0.01, well before the first peak in |C(i)|, the quantum 
surface of section appear essentially localized on the fixed 
phase point corresponding to the horizontal PO, with 
just some probability sticking out of this region along the 
manifolds (also plotted in the figure). For T=0.03 this 
latter effect is more pronounced, and the pattern along 
the linearized part of the stable and unstable manifolds 
is fully developed. At the same time, the corresponding 
low resolution spectral envelopes It(E) (continuous line 
in the bottom tier) define, with increasing precission, the 
clump of eigenstates contributing to this scar ||. Past 
this time, i.e. for T=0.06, well beyond the first (true) re- 
currence of the horizontal PO, the original packet has had 
the opportunity to enter into the truly nonlinear part of 
the periodic trajectory, exploring regions of phase space 
far apart from the fixed point, where the dynamical flux 
gets more complicated. However, we see that the av- 
eraged wave function is well confined on the manifolds, 
even following the sharp kinks that they present. This 
process continues for T=0.1, value at which the local- 
ized wave function appears covering practically all these 
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invariant structures of phase space. In this case the res- 
olution in the spectra is higher, and the corresponding 
spectral envelope define several bands, four in our case. 
These higher resolution bands correspond to the in- 
teraction with other POs of longer periods, as will be 
discussed in the next section. The signatures of these 
longer POs are seen in the configuration space plots. 
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same, if one takes into account the fact that all length 
features scale as fc -1 / 2 . This scaling is the origin of the 
better definition of the focalization effect existing in the 
configuration space wave functions at higher energies. 

Using Gutzwiller trace formulae, similar smoothed 
wave function has recently been calculated |l9|, also 
showing that the invariant hyperbolic classical structure 
is contained in quantum mechanics. In this case no dy- 
namical information was used, with the result that many 
more (of the order of 10000) states were needed in the 
summation in order to get the same kind of localization 
that we are showing in this paper. 

In this section we investigate the dynamics of the lo- 
calized wave functions described in the previous section, 
with the purpose of disentangling the complexity of the 
quantum dynamics in the problem that we are consid- 
ering. As stated in the Introduction, our aim is to pro- 
vide an explanation of the quantum dynamics of highly 
chaotic systems in terms of classical invariant structures. 
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FIG. 3. (Top) Localized (squared) wave functions calcu- 
lated using Eq. (M), (middle) Husimi based quantum Poincare 
surface of section, and (bottom) spectra (7^ as a stick spec- 
trum and It in solid line) for: (a) i?o=400 and T=0.1, and 
(b) £ =40267 and T=0.01. The invariant manifolds of the 
horizontal periodic orbit have also been drawn superimposed 
in the middle pannel. 



IV. DYNAMICS OF THE LOCALIZED WAVE 
FUNCTIONS 

Another aspect of this problem which is worth consid- 
ering is how the process described above is affected when 
the transition to the semiclassical limit is considered. For 
this purpose we have repeated the previous calculations, 
launching the packet from the same point as before but 
with different values of the energy. The results are shown 
in Fig. [|. In it we present the localized wave functions, 
Husimi based quantum Poincare surfaces of section, and 
spectra for £ , =400 and T=0.1 [column (a)], and also for 
£0=40267 and T=0.01 [column (b)]. These values of the 
smoothing times correspond to the periods of the hori- 
zontal PO at the selected energies, so that there are to 
be compared to those in Fig. || (b). As can be seen, the 
structure exhibited by these functions is essentially the 
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FIG. 4. Modulus of the autocorrelation function corre- 
sponding to localized wave functions calculated for: (a) 
#0=400 and T=0.1, and (b) #0=3700 and T=0.03. 

In Fig. U we present the autocorrelation function calcu- 
lated using the (non-stationary) localized wave functions 
[Eq. (0)] of Fig. | (a) and Fig. § (b). The results of part 
(a), corresponding to £o=400 and T=0.1 , show a series 
of recurrences at time intervals approximately equal to 
0.38. In each of them, practically all the initial proba- 
bility is recovered, and the process shows no evidence of 
dispersion up to the times considered in our calculations. 
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t=0.0 t=0.Q42 t=0.QS4 t=0.127 t=0.169 




FIG. 5. Snapshots of the dynamical evolution in phase space of the localized wave function corresponding to Fig. |3| (a). The 
manifolds of the horizontal orbit are plotted in solid line. The linearized manifolds of the orbits (1), (2), (3) and (4) of Fig. ^ 
are included in white lines (in the case that the Husimi function is localized on them). 



This phenomenon has already been discussed in 
Rcf . ji3 , where Gaussian wave packets "stretched" along 
the horizontal PO were used, and necessary conditions 
for the existence of revivals were given. As indicated be- 
fore, in our case we use initial functions that are better 
adapted to the dynamics of the orbit, thus showing an 
improved revival behaviour. In this respect, it is impor- 
tant to note that, opossed to which it was found in p2[ , 
our functions show revivals, in this wavelength regime, 
for an ample range of energy values and also for other 
POs. 

From the quantum mechanical point of view these re- 
currences are easy to interpret, since they are the re- 
sult of the beating frequency originated from the two 
(almost exclusive) eigenstate contribution to the initial 
packet (see spectrum in the bottom part of Fig. ||) . On 
the other hand, at the classical level, the phenomenon is 
more difficult to understand. To further investigate these 
recurrences from this point of view, we show in Fig. || a 
series of snapshots in phase space covering the period of 
time elapsed up to the first recurrence of \ip E °). Notice 
that due to the symmetry of the figure we have only plot- 
ted one half of the complete phase space of the billiard. 
To help interpreting this figure, the position of the fixed 
points and associated manifolds, corresponding to some 
short POs related to the horizontal one, have also been 
included. These orbits are presented in Fig. ^. For the 
sake of clarity, the labelling of the POs was not included 
in Fig. ||, but it can be found in Fig. 0. 

In the first three images of Fig. || we see how the origi- 



nal probability density distribution elongates, stretching 
along the unstable manifolds of the horizontal PO and 
contracting along the stable ones. Afterwards, the pro- 
cess continues, populating the regions around the fixed 
points corresponding to orbits (3), (4) and (1),(2) p0[ , 
so that at t=0. 127-0. 169 a noticeably accumulation of 
probability around these points has taken place. Notice 
also the effect of the flux along the stable manifolds of 
these POs, which at i=0.127 is narrowing the elongated 
shape of the distribution along some specific lines, and at 
t=0.169 has split the packet in several pieces. At t=0.211 
this process reaches its maximum, and the probability re- 
maining in the initial region is practically zero. After this 
point, the evolution continues in a similar way with the 
probability returning gradually to the vicinity of the ini- 
tial point, although obviously this time following a path 
going along the stable manifolds. This second half of the 
dynamics takes place following the same series of inter- 
mediate steps. 

In Fig. U (b) the autocorrelation function for the initial 
state of Fig. @ (b) is shown. In this case the recurrence 
pattern is much more complex, but it is nevertheless no- 
ticeable the existence of a fairly good number of recur- 
rences in which the overlap with the initial distribution 
is greater than 50%. The difference with the example 
presented in the part (a) of the figure is clear. In this 
case, the initial wave function is formed by a substan- 
tially greater number of contributions [see the associated 
spectrum at the bottom of Fig. || (b)], and then the prob- 
ability that all off them get in phase is much smaller. 
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FIG. 6. Some periodic orbits of the desymmetrized stadium 
billiard related to the horizontal one. See text for details. 

The corresponding evolution in phase space is pre- 
sented, up to the first big recurrence at t ~ 0.47, in 
Fig. |^. At this higher energy we see that the mechanism 
for the dispersion of the probability density is practically 
the same, although there exists certain differences worth 
commenting. In the first place, the time scale for the 
dispersion is obviously faster, i.e. compare the extent of 



the evolution in the first snapshot at t=0.059 with the 
third one at t=0.127 of Fig. [| In the second place, the 
probability distribution appears in this case much more 
localized on the participating invariant classical struc- 
tures for all values of time considered. 




FIG. 7. Fixed points and associated manifolds of some pe- 
riodic orbits in Fig. [|The unstable (thick full line) and stable 
(thick dashed line) manifolds of the horizontal periodic orbit 
are also represented. 



t=0.0 t=0.Q59 t=0.110 t=0.162 t=0.213 




FIG. 8. Snapshots of the dynamical evolution in phase space of the localized wave function corresponding to Fig. H (b). 
The manifolds of the horizontal orbit are plotted in solid line. The linearized manifolds of the orbits (5) and (6) of Fig.lq are 
included in white lines at t = 0.213. 
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Finally, at certain values of t, the probability show 
peaks which localize on new POs, not important for the 
dynamical evolution at i?o=400 shown in Fig. |^. For ex- 
ample, at t — 0.213 the probability density localize on or- 
bits (5) and (6) of Fig. ||. This effect is not unexpected, 
since it is reasonable to admit that as we increase the 
energy, going towards the semiclassical limit, more POs 
are necessary to explain the quantum dynamics of our 
system, as discussed in Ref. ||. 



V. SUMMARY AND CONCLUSIONS 



quantum mechanics of very chaotic Hamiltonian systems 
in terms of simple classical structures. 
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The introduction of Gutzwiller's trace formulae g] in 
1971 provided a method to understand the quantum 
mechanics of classically chaotic Hamiltonian systems in 
terms of POs. Its main drawback is the exponential 
growth in the number of orbits involved as the Heisen- 
berg time, Th, increases, and many efforts have been 
devoted to the development of resummation techniques 
that improve its convergence. 

Our work considers an alternative approach, initiated 
with Refs. based on the use of special functions well 
localized on a few POs and its manifolds, the number of 
which grows only linearly with Th- To construct these 
localized functions we have used the method described in 
Rcf . H , based on the dynamics of a coherent wave packet 
averaged along a given PO and its neighborhood. 

In this paper we have studied in detail some important 
characteristics and properties of these non-stationary 
scarred functions in the case of the horizontal PO in the 
desymmetrized Bunimovich stadium billiard, a paradig- 
matic example of chaos. In particular, we have consid- 
ered how the localization properties, both in configura- 
tion and phase space, of these functions calculated at 
energies, Eq, close to the Bohr-Sommerfeld quantization 
conditions vary with the averaging time, T. Our results 
show that the probability density in phase space extends 
along the manifolds of the scarring PO, even when the 
associated dynamics have left the linear regime. In the 
second place, the transition to the semiclassical limit has 
also been investigated by simultaneously changing E 
and T. As it takes place, all length features get better 
defined, for example the focal points corresponding to the 
hyperbolic structure around the fixed point. Finally, we 
have also studied the dynamics of these non-stationary 
wave functions. When the corresponding probability 
density in phase space is followed a dispersion along the 
manifolds is observed. In this process a building up of 
probability peaks takes place on a few POs related to the 
original one (the horizontal PO in our case). The recur- 
rences in the corresponding autocorrelation function can 
be very important, and we have even shown one example 
in which they are particularly pronounced. 

To conclude, we think that this work provides valuable 
information that can help in the disentanglement of the 
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